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We probe the flow of two dimensional foams, consisting of a monolayer of bubbles sandwiched 
between a liquid bath and glass plate, as a function of driving rate, packing fraction and degree 
of disorder. First, we find that bidisperse, disordered foams exhibit strongly rate dependent and 
inhomogeneous (shear banded) velocity profiles, while monodisperse, ordered foams are also shear 
banded, but essentially rate independent. Second, we introduce a simple model based on balancing 
the averaged drag forces between the bubbles and the top plate, Fbw and the averaged bubble- 
bubble drag forces Fbb, and assume that Fbw ~ v 2 ^ and Fbb ~ (dyv) 13 , where v and (d y v) denote 
average bubble velocities and gradients. This model captures the observed rate dependent flows for 
P ~ 0.36, and the rate independent flows for j3 « 0.67. Third, we perform independent rheological 
measurements of Fbw and Fbb, both for ordered and disordered systems, and find these to be fully 
consistent with the forms assumed in the simple model. Disorder thus modifies the exponent /3. 
Fourth, we vary the packing fraction (f> of the foam over a substantial range, and find that the flow 
profiles become increasingly shear banded when the foam is made wetter. Surprisingly, our model 
describes flow profiles and rate dependence over the whole range of packing fractions with the same 
power law exponents — only a dimensionless number k which measures the ratio of the pre-factors 
of the viscous drag laws is seen to vary with packing fraction. We find that k ~ ((f) — <^ c ) _1 , where 
4>c ~ 0.84, corresponding to the 2d jamming density, and suggest that this scaling follows from 
the geometry of the deformed facets between bubbles in contact. Overall, our work shows that the 
presence of disorder qualitatively changes the effective bubble-bubble drag forces, and suggests a 
route to rationalize aspects of the ubiquitous Herschel-Bulkley (power law) rheology observed in a 
wide range of disordered materials. 

PACS numbers: 47.15.gm, 47.57.Bc, 83.50.Lh 



I. INTRODUCTION 

Foams, which are dispersions of densely packed gas 
bubbles in a liquid, exhibit an intricate mix of elastic, 
plastic and viscous behavior reminiscent of the mechan- 
ics of other disordered materials such as colloidal sus- 
pensions, granular media and emulsions [TJ [2j [3j HJ [6] . 
When left unperturbed, foams jam into a meta-stable 
state where surface tension provides the restoring force 
underlying their elastic response for small strains [TJ[IJ[7]. 
Under continuous driving the foam starts to flow, and the 
viscous dissipation that arises in the thin fluid films that 
surround the gas bubbles becomes important. Macro- 
scopically, the steady state rheology of foams exhibits 
shear thinning, and the stress r as function of strain rate 
7 is generally non- linear, often taking a Herschel-Bulkley 
form: r = Ty + C17' 3 , where Ty denotes the yield stress, 
and where the viscous stress ry = r — Ty scales nontriv- 
ially with the strain rate -y pO [HI HI El EH1 HU H21 H31 CU ■ 
In addition, in many situations, the flow is inhomoge- 
neous and localizes in a shear band [TJ [T2 [HJ [T5J [TB] . 

In an earlier Letter [T3] we experimentally probed the 
flow of disordered, bidisperse 2d foams which are trapped 
between the fluid phase and a top-plate. The 2d nature 
allows for direct imaging of the bubble dynamics and of 
the shear banded flow profiles in this system. Combin- 
ing measurements of the flow profiles with rheological 
measurements, we established that viscous interactions 
between neighboring bubbles scale differently with ve- 
locity gradients than the effective viscous interactions at 



the global scale. We captured the rate-dependent shear 
banding exhibited by our system in a nonlinear drag force 
balance model. Here we expand on these findings, dis- 
cuss new results for the effect of varying the wetness of 
the foam, and provide extensive additional evidence to 
support our main conclusions. 

To understand the rheology and shear band formation 
in our system, three ingredients need to be described and 
combined appropriately: (i) Interactions with the top 
plate, (ii) Local bubble interactions. (Hi) Disorder. 

Top plate — In recent years, a variety of studies have 
addressed the formation of shear bands in (quasi) two- 
dimensional foams, consisting of a single layer of macro- 
scopic bubbles. Such single layers can be made by freely 
floating the bubbles on the surface of a surfactant so- 
lution ("bubble raft") [H EJ JTB], by trapping them 
between a top glass plate and the surfactant solution 
("liquid-glass") [3 HHJ QU [20j [21] , or by trapping them 
between two parallel glass plates (Hele-Shaw cell) [To] . 

In a seminal paper by Debregeas et al. [TS], a bidis- 
perse foam in Hele-Shaw Couette cell was sheared and 
narrow shear banded flow profiles where obtained |15j . 
While initially it was believed that for slow flows, the 
effect of the viscous drag forces exerted by the confin- 
ing glass plates would be negligible [TJ [22], these drag 
forces have turned out to be crucial. First, Couette ex- 
periments in bubble rafts found completely smooth flow 
behavior |12j . Second, in experiments where a monodis- 
perse foam was linearly sheared with and without confin- 
ing glass plate on top [TS], one observes smooth velocity 
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profiles for the bubble raft but highly shear banded flows 
for the liquid-glass geometry. The precise connection be- 
tween the drag forces due to the confining plates and the 
occurrence of shear banding in confined foams is still a 
subject of debate [2"2l |2"3] . 

A simple continuum model that balances the top plate 
drags and the inter-bubble drags (modeled with a Bing- 
ham constitutive relation) captures both the rate inde- 
pendence and exponentially localized shear bands seen in 
the linear liquid-glass cell [TU [H] — here we will build on 
and extend this model to capture the experimentally ob- 
served nonlinear, rate dependent rheology of disordered 
foams. 

Local interactions — At the microscopic level, bubble 
interactions are a combination of elastic repulsion, typi- 
cally harmonic for small deformations ES, 26 , 27 , 28], and 
nonlinear viscous drag forces [T31 [Til ®§\ |3tJl [5Tj |32l 133] . 
Such drag forces arise when two bubbles slide past each 
other or when a bubble slides past a solid boundary. The 
viscous drag forces originate in the thin films that sur- 
round foam bubbles, and have recently received renewed 
attention [T31 [221 EHJ EH EH]- Already for a single 
bubble sliding past a solid wall, Bretherton showed that 
the drag force scales nonlinearly with the bubble velocity 
[T31 [521 134"] , and by analogy one would expect the drag 
forces arising between sliding bubbles to be nonlinear also 
— indeed Denkov et al. recently suggested that a similar 
scaling applies to the viscous drag force between bubbles 

031. 

Here we measure these drag forces directly by rheolog- 
ical experiments where two rows of ordered bubbles are 
sheared past each other. 

Disorder — Foam flows are disordered and intermit- 
tent at the multi-bubble scale [5j [T3J HH EH [33 [36] . 
For such disordered systems, the affine approach, where 
one simply scales up local elastic or viscous interactions, 
often fails to describe the macroscopic behavior — this is 
by now well established for shear deformations in gran- 
ular and foam- like systems [371 13H1 EH], and a similar 
picture is emerging from simulations of the flow of vis- 
cous particles [351 HOI EH 113 US]- In the present work 
we present strong experimental evidence for the failure 
of the affine approach to describe drag forces in flowing 
systems. 



Outline 

In this paper, we describe an experiment in which we 
have linearly sheared a 2d foam and we disentangle the 
roles of the top plate, the local bubble interactions and 
the disorder, as well as the role of the wetness of the 
foam. 

In section ITT] we describe our experimental setup. In 



hibit rate dependent flow profiles, which become increas- 
ingly shear-banded for large 7 a . We capture our find- 
ings in a model in which the time-averaged drag forces 
between bubble and top plate, F{, w , and between neigh- 
boring bubbles, are balanced. While the continuum 
limit of our model is similar in spirit to [24] . the cru- 
cial new ingredient is nonlinear scaling laws for the wall 
drag and the bulk stress — these nonlinear scalings are 
essential for capturing the observed rate dependence. 

In section |IV] we probe the scaling of the drag forces 
by independent rheological measurements, allowing us to 
directly probe the role of disorder by comparing the rhe- 
ology of small ordered and larger disordered bubble rafts. 
We find the averaged drag forces in the disordered foam 
to scale differently from the local drag forces between 
individual bubbles, which we have measured at high res- 
olution and analyze in a novel way. 

In contrast, for monodisperse, ordered foams the lo- 
cal, averaged and top-plate drag forces all scale similarly, 
causing rate-independent flows similar to those seen by 
Wang et al. [18], and we discuss these in section [v] 

In section [yT] we further probe the connection between 
the viscous drag forces at the bubble scale and the bulk 
viscous forces by performing additional linear shear ex- 
periments over a range of packing (air) fractions <f>. We 
find that the contribution of averaged bubble-bubble drag 
forces vanishes algebraically as <f> — <j) c , when the packing 
fraction is decreased towards a critical value <f> c , which we 
identify with the (un)jamming density — <f> c ss 0.84. We 
relate the vanishing of the averaged bubble-bubble drag 
forces at (j) c to the vanishing overlap between bubbles at 
unjamming. 

The simple clastic interaction (typically harmonic for 
small deformations E3 EH1 EH EE]) and the absence of 
solid friction make static packings of foam bubbles em- 
inently suited to compare to simulations of the popular 
soft frictionless sphere model [38] HH [35] . Our work illus- 
trates the great potential of foams to elucidate the flow 
behavior of simple systems near jamming [40 ] liTl |42"1 143] . 



section III we present experimental results for flow pro- 
files for a range of strain rates and span- wise widths of our 
system. We find that the flow depends crucially on the 
applied strain rate j a '- disordered, bidisperse foams ex- 



II. EXPERIMENTAL DETAILS 

In this section we describe in detail a novel experimen- 
tal setup to induce linear shear flow in two-dimensional 
foams. We also detail the analysis techniques used to ex- 
tract velocity profiles, and discuss measurements which 
show that coarsening and fluid drag can be neglected. 



A. Setup 

We create foam bubbles on the surface of a reservoir 
of soapy solution (of depth 3.5 cm), consisting of 80% by 
volume dcmincralized water, 15% glycerol and 5% Blue 
Dawn dishwashing agent (Proctor & Gamble), by bub- 
bling nitrogen through the solution via syringe needles of 
variable aperture. We measure the bath surface tension 
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a with the pendant drop method [IB] and find a = 28 
± 1 mN/m. We measure the dynamic viscosity rj with 
a Cannon Ubbelohde viscometer and find 77 = 1.8± 0.1 
mPa.s. 

Fig. 1 shows our experimental setup: the bubbles 
are contained inside an aluminum frame (400x230 mm) 
which is leveled with the liquid surface and which sup- 
ports glass top plates to which the bubbles bridge once 
they are in place. The top plates consist of three adja- 
cent glass plates with slits to accommodate two PMMA 
wheels of radius 195 mm and thickness 9.5 mm which 
drive the flow. The vertical gap between the liquid sur- 
face and the glass plates can be varied to control the 
packing fraction of bubbles <fi. 

The wheels, which are grooved to provide a no-slip 
boundary for the bubbles, can be lowered into and raised 
out of the bath through the slits. The wheels are con- 
nected to two Lin Engineering stepper motors, each 
driven by micro-stepping driver, and are rotated in op- 
posite directions. At any point along the line where the 
wheels contact the foam bubbles the horizontal compo- 
nent of the driving velocity is a constant (see Fig. lb). 

We obtain our data from the central 60 mm of the 
shearing region — marked by the horizontal lines in 
Fig. [TJa) — to avoid effects caused by the recirculation 
of the foam at the edges of the wheels. In this central 
part no motion is observed due to the vertical component 
of the radial velocity. At the edges of the slits, bubbles 
do leave the system, while being pinned to the wheels. 
This does not result in holes in the foam layer, either 
because at high driving velocities the bubbles reenter the 
system before rupturing while traveling on the wheel, or 
because at low velocities bubbles from outside the shear- 
ing region are pushed inwards due to the bubble surplus 
at the edges. The resulting driving velocity gives rise to 
a global strain rate 7 = 2vq/W, where W denotes the 
gap between the wheels, which we vary between 5 and 10 
cm. 



B. Imaging and Analysis 





FIG. 1: (Color online) (a) Schematic topview of the experi- 
mental setup. W represents the gap width and the two hor- 
izontal lines indicate the edge of the region over which the 
velocity profiles are calculated. The red curve depicts one 
such profile, (b) Sideview of the shearing wheels. The slits in 
the glass plate are drawn for clarity. That the in-plane com- 
ponent of the motion of the boundary is constant can be seen 
as follows: by trivial geometry, we obtain that vo=u)ri cos</>, 



but since n = 



COS <f) ' 



at any point along the along the contact 



line of 230 mm, the layer of bubbles is sheared with a driving 



velocity vo — lo- 



ujr . 



(c) Experimental image of 



part of the foam, the scalebar represents 5 mm. 



in the bubble displacements, we take 1000 frames per 
run, corresponding to a strain of 4 for a 5 cm gap, as we 
are interested in averaged velocity profiles. We pre-shear 
the system before taking data, so that a steady state is 
reached. 

We obtain the velocity profiles both through particle 
tracking and a Particle Image Velocimetry-like technique, 
where for each y-value, we calculate the cross-correlation 
(C n ) 2 between the corresponding image line in the P n {x) 
of length m and the same image line P n+ i{x) in the next 
frame shifted by an amount r: 



We wish to characterize the average flow in the re- 
direction as a function of the span-wise coordinate y. 
The average velocity profiles are obtained from a series 
of images which we record with an 8 bit Foculus BW 432 
CCD camera (1280x1024 pixels) equipped with a Tam- 
ron 28-300 telezoom objective. In the images, 1 pixel 
corresponds to approximately 0.1 mm. To optimize the 
brightness and obtain images in which the bubbles ap- 
pear as circles, the foam is lit laterally by two fluorescent 
tubes, each driven by high frequency ballasts to prevent 
flickering in the images. The bottom of the reservoir is 
covered with a black plate to improve contrast. Typical 
images are shown in Fig. [2] 

The frame rate is fixed such that the displacement at 
the wheels is fixed at 0.15 mm between frames. Since 
the flow is strongly intermittent, with large fluctuations 



(C„(r)) 2 = V P n (i)P n+ i(i + T). 



=0 



(1) 



We can then proceed in two ways. The first method is 
to add up all cross-correlations from all frames for each 
y-value, and calculate the average displacement Aa; per 
frame by fitting a parabola p n {r) to the resulting sum 
of cross-correlations and taking the peak value of that 
parabola: 



'999 n 

Aa;(y) = max I VV>(r) 



(2) 



In the second method we fit a parabola to each cross- 
correlation separately and obtain the average displace- 
ment by averaging the maxima of all individual parabo- 
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FIG. 2: (Color online) Images of sheared regions for both (a) 
monodisperse and (b) bidisperse foams. Shear is indicated by 
the arrows. The highlighted area is where data analysis is 
performed on. 



las: 



Ax(y) = (max(p„(r))). 



(3) 



By comparing to average velocity profiles obtained by 
particle tracking [33], we find that the latter procedure 
gives the closest match to the tracking velocity profiles, 
and we have employed this procedure throughout. We 
thus obtain both spatially (in the x-direction) and tem- 
porally averaged velocity profiles. Despite the intermit- 
tent character of the flow, we obtain smooth reproducible 
velocity profiles. 



I • t= hn 
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FIG. 3: (a) Size distribution and coarsening over the duration 
of an experimental run for bidisperse foams, (b) Flow at 
the liquid surface in the absence of bubbles, as imaged by 
depositing silver powder. Inset: same profile on lin-log scale, 
showing exponential decay away from the boundaries. 



We furthermore measure the velocity profile of the liq- 
uid surface itself at the same fluid level as in the foam 
experiments by imaging the flow of silver particles that 
were sprinkled on the liquid surface, see Fig. [SJd. We 
observe a steeply decreasing velocity profile at the fluid 
surface, which implies that even if the fluid drag were 
of the order of the other drags acting on the bubbles, it 
would not significantly alter the flow profiles except near 
the wheels. 

We thus conclude that the bubble size distribution is 
essentially constant during an experimental time frame, 
and that the dominant drag forces are those between bub- 
bles and top plate, and those between contacting bubbles. 



III. LINEAR SHEAR OF TWO DIMENSIONAL 
FOAMS 

In this section we explore the rate dependent shear 
flows in our system experimentally. By fitting our ex- 
perimental data to a nonlinear drag force balance model, 
we deduce the dependence of the averaged bubble-bubble 
and bubble- wall drag forces as function of the local strain 
rate and velocity. 



C. Coarsening and fluid drag 

To characterize the amount of coarsening we measure 
the bubble size distribution by measuring the surface area 
of the bubbles in the images. We obtain well defined 
size distributions which show little coarsening over the 
duration of the runs, which corresponds to about 2 hours 
(Fig. 3a). 

We have checked that the drag on the foam bubbles 
due to flow of the bulk liquid underneath is negligible 
by measuring the velocity profile of bubbles floating on 
a very shallow layer of bulk fluid. In this case the fluid 
surface velocity is decreased due to the no-slip boundary 
condition at the reservoirs' bottom. We do not, within 
experimental uncertainty, observe a change in the exper- 
imental velocity profiles in this geometry. 



A. Flow of disordered foams 

We measure averaged velocity profiles in disordered 
twodimensional foams. These foams are produced by 
bubbling a fixed flow rate of nitrogen through syringe 
needles of 2 different inner diameters, such that bubbles 
of 1.8 ±0.1 and 2.7 ±0.2 mm result (at 59-41 number 
ratio). The bubbles are gently mixed with a spoon un- 
til a disordered monolayer results. For gap widths of 5, 
7 and 9 cm, we drive the foam at 6 different velocities, 
spanning 2.5 decades: v = 0.026, 0.083, 0.26, 0.83, 2.6 
and 8.3 mm/s. 

Note that we perform the sweep in driving velocities 
from fast to slow and that we pre-shear the system for 
one full wheel rotation, to start with bubbles covering 
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FIG. 4: (Color online) (a-c) Flow profiles for a gap width W = 5 (a) 7 (b) and 9(c) cm. From black to light grey, vq = 0.026 
mm/s, 0.083 mm/s, 0.26 mm/s, 0.83 mm/s, 2.6 mm/s and 8.3 mm/s. For all gap widths we observe that the localization near 
the driving wheels increases for increasing driving velocity. For clarity the profiles are each offset vertically by 0.5xv/vo. Solid 
red curves: fits to the drag force balance model of section III B. (d) Profiles at 2.6 mm/s for all three gap widths. Regardless 
of the gap width all profiles decay at the same rate, (e) Examples of profiles and fits on a log-log plot, highlighting the linear 
tails of the profiles, vo — 0.026, 0.26, 2.6 mm/s, W = 5 cm. (f) As explained in section III C 1, the minimum in /3 is found by 
calculating the variance ^/(Sk 2 )/k 2 of k over all 18 profiles depicted in (a-c). The minimum in the variance is seen at /3 = 0.36 
— see section III C. 



the wheel [J7| and ensure that we have reached a steady 
state. To fix the packing fraction, we fix the gap between 
glass plate and liquid surface at 2.25 ± 0.01 mm. We have 
measured (see section VI) that for this gap the packing 
fraction is <j) = 0.965 ± 0.005. 

Results are plotted in Fig. |3Ja-c) : the profiles exhibit 
shear banding, and for all gap widths the profiles become 
increasingly shear banded at increasing driving velocities. 
The slowest runs at W — 5 cm yield essentially linear ve- 
locity profiles. We suggest that these shapes are due to 
the small gap width, which results in overlapping shear 



banded profiles resembling a linear profile, and in what 
follows, we will present a model that supports this con- 
clusion. 

In Fig. Qd) we plot velocity profiles for a driving ve- 
locity of 0.26 mm/s for all three gap widths together, 
which clearly show that for all widths, the velocity pro- 
files decay similarly. Fig. Qd) thus suggests that in this 
experiment the driving velocity at the edges, instead of 
the overall shear, sets the velocity profiles, and that the 
local response to forcing will provide the key towards un- 
derstanding the shape of these profiles. Note finally that 
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the profiles do not exhibit significant slip with respect to 
shearing wheels, except for the fastest runs, where the 
slip is less than 20 %. 



B. Model 

We now propose a model to account for the shear band- 
ing behavior discussed above, by considering the balance 
of the averaged viscous drag forces. 



1. Drag forces on individual bubbles 

The drag force on a single bubble that slides past a 
solid wall was first investigated by Bretherton [53] and 
has recently received renewed attention |13j |29j [30j [3TJ 
1521 H5] . The crucial finding is that F bw , the drag force 
per bubble sliding past a solid wall, scales as 



Fbw = /ta(Ca) 2 / 3 = f bw (rjv/a) 



2/3 



(4) 



with r\ the bulk viscosity, a the surface tension, a 
constant with dimensions of force and Ca the capillary 
number. Typically fbw °c or c , with r c the radius of the 
deformed contact between bubble and wall [48]. For bub- 
bles in a soapy solution, the 2/3 scaling with Ca only 
holds for surfactants that are mobile [13]. Results from 
[4"9"] strongly indicate that this is indeed the case for our 
surfactant Dawn, as we will confirm below. 

The drag force between 2 bubbles sliding past each 
other, Fbb, has not received much attention up to now, 
although [50] provides indirect evidence that it scales like 
Fbb °f (r/ACa)^ , with ACa = r)Av/a. In a very recent 
Letter it is explicitly shown that, for ordered bubble mo- 
tion Fbb scales indeed as (ACa)'' [33]. The authors find 
(, = 0.5, although various physico-chemical peculiarities, 
as well as the range of Ca one measures in, can alter this 
exponent. 

Taking all of this into consideration, it seems reason- 
able to assume that: 



F bb = fbb (rjAv/a) 



C 



(5) 



While the dissipation leading to Fb w occurs at the 
perimeter of the flattened facet [T3] — hence the pref- 
actor fbw oc crr c — fbb scales oc <tk 2 , with k c the radius 
of the deformed contact between bubbles, thus reflecting 
the different physical mechanism behind this scaling 33J . 



2. Stress balance 

We divide our shearing region in lanes labeled i and 
assume that on every lane the time-averaged top plate 
drag per bubble F % bw balances with the time-averaged 
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FIG. 5: (Color online) Illustration of drag balance model. The 
shear region is divided in lanes labeled i which all experience 
drag forces due to the top plate and due to both neighboring 
lanes, (b-c) Illustration of the films around which the viscous 
drag forces act. 



viscous drag per bubble due to the lane to the left (F bb ) 
and right (F bb ), see Fig. 



pi+i 

r bb 



F 



bw 



F 



bb 



0. 



(6) 



We assume that the averaged drag forces scale similarly 
to the local drag forces. For the averaged bubble-wall 
forces we assume: 



Fl 



fUw*/*) 2 ' 3 , 



(7) 



while for the averaged bubble-bubble drag forces we as- 



sume: 



Fl = fY + faiivMiv'-v'- 1 )] 1 
Fit 1 = fY + fbbiivHiv^-v*)] 







(8) 
(9) 



Here f bw and fbb are material parameters with dimension 
of force, which will be measured by rheometry in section 
IV below. Finally, fy represents a yield force in the inter- 
bubble drag, to remain consistent with rheometrical data 
presented later on and to reflect the elastic barrier bub- 
bles have to overcome before they slide past each other. 
Note that the velocities v % denote the averaged velocities 
in the x-direction — the crucial assumption is that the 
relation between the averaged drag forces and the aver- 
aged velocities is simple and can be expressed by a single 
power law. 

We do not know if and how the conjectured forms for 
the averaged forces can be derived from the non-averaged 
forces Eqs. Q and (|5| since due to the intermittent and 
disordered bubble motion, the instantaneous bubble ve- 
locities are fluctuating and not necessarily pointing in 
the x-direction. For example, there is no a priori rea- 
son for the exponents £ and to be equal and in fact 
our data strongly indicate that they are not. The best 
justification for Eqs. (7|9l is a posteriori — the result- 
ing model describes the data well. Note that the bars in 
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Eqs. (7|9l express an average over disorder, in the sense 
that these quantities are measured in highly disordered, 
intermittent flows. 

Inserting the expressions from Eqs. (7|9l into Eq. ^ 
and defining k — fb w / fbb we arrive at: 



(10) 

Note that the yield drag contributions fy cancel, which 
is a particular advantage of the linear geometry we work 
in. The model predicts flow profiles for arbitrary width 
and driving rate, once the parameters (3 and k are fixed. 



Fits 



1 . Procedure 

We compare all 18 runs to solutions of the model. We 
focus on the central part of the data where < 3/4i>o 
to avoid the edge effects near the shearing wheels (for in- 
stance the bumps in the low-velocity profiles in Fig. Qa) 
and the slip with respect to the wheel in the fast runs). 
We numerically integrate Eq. (10 1 from y — 0, where 



D. Continuum Limit 



v = 0, to the y value for which v — 3/4 • vq, while keeping 
(3 and k fixed. The drag force balance should govern the 
shape of the velocity profiles for all driving rates and gap 
widths. Therefore we determine for fixed /3 the k values 
that fit the flow profiles best. The k values exhibit a 
systematic variation that depends on the value of (3. We 
quantify this variation by computing the relative variance 
■>J {8k 2 ) /k 2 and by repeating the procedure for a range 
of [3, we obtain a plot of the variance as a function of /?, 
see Fig. Qf). From this graph, we determine the value 
for which the variance is minimized as (3 = 0.36 ± 0.05. 



2. Results 

Fixing now k — 3.75 and [3 = 0.36, we capture the 
shape of all data sets with high accuracy . The resulting 
model profiles are plotted in Fig. |4ja-c) , and we see that 
for these values all velocity profiles are adequately fitted 
except for the slowest runs at W = 5 cm. We attribute 
these deviations for small W to the observation that edge 
effects extend further into the shearing region for small 
gaps. 

Note that the model profiles exhibit linear tails, see 
Fig.[4je), and that the experimental velocity profiles ex- 
hibit approximately the same behavior. We conclude 
that both the experimental and model profiles do not 
decay exponentially, in contrast with results found in pre- 
vious studies H21 HH] ■ 



The continuum limit of Eq. ( 10 1 can be written as: 

2/3 



dr 
dy ' 



(11) 



r = m /,,,, ( UMl2]\ i . (12) 



Hence, the top plate drag can be considered as a body 
force and the inter-bubble drag force as the divergence 
of a shear stress r, where Ty is an undetermined yield 
stress. Eq. [12] is the constitutive equation for a Herschel- 
Bulkley fluid [51 , and we can now associate the averaged 
bubble drag force scaling at the local level with the power 
law scaling of the viscous stress in the Herschcl-Bulklcy 
model. 

Note that (3 — 0.36 is similar to the power law index 
n = 0.40 found for the bulk rheology of three-dimensional 
mobile foams [31 [T3] and to the values n = 0.33 and 
n = 0.45 found for two dimensional bubble rafts in a 
Taylor-Couette geometry in [IB"] . 

The fact that the yield stress does not play a role for 
our velocity profiles can now be understood in two ways: 
at the continuum level, since it is a constant it vanishes 
after taking the divergence of the shear stress, at the bub- 
ble level, even though we include a yield force in Eqs. 0- 
|8j, the contributions from both neighboring lanes cancel 
in Eq. ( 10 1. Finally, notice that the continuum equations 



can easily be solved in terms of hypergeometric functions 
[52]. 



IV. RHEOMETRICAL DETERMINATION OF 
VISCOUS FORCES IN TWO-DIMENSIONAL 
FOAMS 

In this section we will investigate the viscous forces 
that act at the bubble scale by rheometry, to test and val- 
idate the assumptions for the scaling of the bubble-wall 
drag and the viscous friction inside the foam expressed 
in Eqs. ppf . We use an Anton Paar DSR 301 stress 
controlled rheometer, which can also operate in strain 
controlled mode. We use the rheometer in strain con- 
trolled mode to investigate Fb w - Moreover, we compare 

measurements, which reflect the actual drag force at the 
single bubble level (Fbb), with measurements of the aver- 
aged viscous drag force on a bubble in a disordered flow 
of foam (Fbb)- 



A. Bubble-wall drag 

We directly measure the bubble-wall friction for foam 
bubbles produced from the soap solution presented 
above, with a method that was introduced in [T3J. We 
load a monolayer of bubbles between two PMMA plates 
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FIG. 6: (Color online) (a) Drag force per bubble exerted on smooth rotated plate as a function of Ca, probed by the total drag 
force of a pinned layer of bubbles on a rotating top plate. The solid line represents Fh w = 0.0015 ±0.0001 • (r/v/a) 2 ^ 3 . The upper 
inset shows a close-up photograph of the rheometrical tool used to measure the bubble-wall drag: the reflection of the flattened 
facets of radius r c used to extract Rq can be seen clearly. The lower inset shows a side view of the experimental geometry, (b) 
Raw torque for ordered and commensurate lanes of bubbles sliding past each other. Notice the huge fluctuations with respect 
to the mean indicated by the horizontal line . The average of the raw data corresponds to the data point in (c) indicated by 
the arrow, (c) Torque averaged over an integer number of rearrangements as a function of ACa for the commensurate case 
(40 bubbles on inner wheel, 40 bubbles on outer wheel)(B), incommensurate case 41/40 (T), incommensurate case 44/40 (a). 
Dashed lines indicate £ = 0.67 resp. £ = 0.75. Upper inset shows a schematic picture of the rheometrical geometry, lower inset 
shows a histogram of the extracted values of the exponent £. The width of the bin indicates the error in f . 



of radius Rp — 2 cm. The bubbles are pinned to the lower 
plate by means of a hexagonal pattern of indentations of 
size O(d), and can slip with respect to the smooth upper 
plate which is connected to the rheometer head, see lower 
inset of Fig. [6ja) . We measure the torque T exerted by 
the bubbles as a function of the angular velocity oj of the 
smooth plate. 

We convert T(uS) to Fb w (Ca) in the following way: 
each bubble exerts a wall stress t w = F^/itRq on 

the smooth plate. We integrate the contribution to the 
torque of this wall stress over the plate: 



T 



T w r2-Krdr 



2r 2 dr. (13) 



If we now assume that Fb w °c [Ca] Q = [^^-j", we 
can immediately read of from the data that a = 0.67, 
see Fig. |6ja), so inserting this expression in the integral 
Eq. (pi) yields: 



T 



2Fh w R, 



3.67 



3.67RI 



(14) 



Since the bubbles are flattened during the measurement, 
we can only measure Rq through the flattened facet r c 
by looking at the reflection of the deformed facet, see the 
upper inset of Fig. ^&). We find r c — 1.59 ± 0.05 mm. 
As the bubble radius is smaller than k _1 we can express 

R Q in terms of r c through Rq = \J\tcK~ 1 [15] . Note 
that this derivation of r c in terms of Rq hinges on the as- 
sumption that the bubbles are not too deformed, which 
is not obvious in the rheometrical geometry, but for lack 



of a more precise relation we use it. We finally rescalc the 
horizontal axis by multiplying lj with r\R p j a . The result- 
ing curve is plotted in Fig. [6]ja): over our measurement 
range (more than three decades) F^ w cx [Ca] 2//3 . 



B. Bubble-bubble drag 

1. Drag at the bubble scale 

To measure the power law scaling of the inter-bubble 
drag we measure the torque exerted by a foam driven 
at a strain rate 7 in a cylindrical Couette geometry, 
which consists of an inner driving wheel, connected to the 
rheometer head, rotating inside an outer ring. The rheo- 
metrical experiments are performed with bubble rafts, 
i.e. foams that are not confined by a top plate, as the 
additional stresses due to the wall would disturb a clean 
rheological measurement. 

Both boundaries are grooved to ensure a no slip bound- 
ary for the bubbles, of which a monolayer floats in the 
shearing region. We start with measuring F^ for the 

ordered case by keeping the gap between the cylinders 
such that exactly two layers of bubbles fit in, see the up- 
per inset of Fig. [6jc) . The inner radius (r,) is 2.5 cm 
and the outer radius (r G ) is 3.0 cm. We deposit bubbles 
of 2.2 mm diameter in the grooves, make sure that all 
bubbles are strictly pinned and remain in their groove, 
and vary the rotation rate lo of the inner cylinder over 
3 decades while measuring the torque averaged over an 
integer number of rearrangement events, see Fig. |6jc). 
We multiply u> by rjri/a to rescale the dimensionless 
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FIG. 7: (Color online) (a) Torque exerted on the inner wheel by a monodisperse foam in a Taylor-Couette geometry with 
Ti = 25mm, r = 70mm, for different bubble sizes as indicated. Fits are to a Herschel-Bulkley model, and power law indices /3 
from these fits are shown in the graph as well. Surprisingly, the yield stress increases with increasing bubble size — see text, 
(b) Same data as in (a) with the yield torque from the fit subtracted. The solid line is a power law with index 0.4. (c) Averaged 
drag force per bubble in a bidisperse, disordered foam. The foam is sheared in a Couette cell of inner radius 5 cm, outer radius 
7 cm (hence a gap of 9 bubble diameters) without a top plate, see inset. We obtain Fbb = fy + /6t(ACa) /3 , with the yield 
threshold f Y « 2.2 ± 0.5 x I0~ 6 N, f bb « 2.5 ± 0.9 x f 0~ 4 N and /3 = 0.40 ± 0.02 (solid line). Open circles are the same data 
with the yield torque obtained from the fit subtracted, which are well fit by a pure power-law with exponent 0.4 (dashed line). 



velocity difference and we divide the torque by T\ and the 
number of bubbles pinned at the inner wheel (e.g. 40) to 
obtain the averaged bubble-bubble drag force per bubble 
in the ordered case. 

We use three different inner wheels: one with 40 
grooves, a second with 41 grooves and a third with 44 
grooves. Since the number of grooves in the outer ring 
is fixed at 40, this allows us to investigate the differences 
between commensurate and incommensurate numbers of 
bubbles in the grooves. 

For the commensurate case, the result is plotted in 
Fig.^b): All bubbles rearrange simultaneously and thus 
the signal reflects the torque exerted on a single bub- 
ble, amplified by a factor of 40. The elastic barrier that 
bubbles have to cross before rearranging is clearly visible 
in the signal. As a result, the torque oscillates tremen- 
dously. Nevertheless, the force per bubble averaged over 
many such events scales with the dimensionless velocity 
difference as a power law with index 0.7, see Fig. [6] (c). 
This value is remarkably close to the exponent found for 
the bubble-wall drag. For these ordered lanes, no signs of 
a yield plateau are observed in the time averaged signal, 
and we believe this is due to the fact that all elastic en- 
ergy that is stored in the bubble deformation is released 
after yielding, so that one measures purely the viscous 
drag. 

For the incommensurate runs, the raw signal looks 
more complex, as rearrangements do not occur simulta- 
neously for all 41 or 44 bubbles. The resulting power-law 
exponents for the averaged drag forces are, however, close 
to the one observed for the commensurate case. In fact, if 
we repeat the measurements for both commensurate and 
incommensurate bubble numbers a multitude of times 
and fit Caf' to the averaged Fw, we find a distribution 



of ^-values around ( = 0.73, see lower inset of Fig. ^c). 
The binsize is similar to the errorbar on each individual 
measurement. 



2. From local to bulk viscous drag 

We observe that the scaling exponent for the viscous 
drag at the bubble scale, £, differs markedly from the 
scaling exponent j3 of the drag forces inside the bulk foam 
as extracted from the velocity profiles, e.g., £ » 0.70 vs. 
/3 = 0.36. We hypothesize this is due to the disordered 
flow in the foam and will provide rheological evidence in 
what follows. 

To perform rheological measurements of the drag 
forces, we employ a Couette cell which has an outer ring 
of radius ro — 7 cm, such that more layers of bubbles can 
fit inside the cell. We first will perform measurements on 
disordered packings of monodisperse bubbles of three dif- 
ferent sizes (1, 3 or 5 mm). We observe that the foam de- 
viates substantially from hexagonal packing during flow 
because the inner radius Ti = 2.5 cm is small, and the 
curvature is large. We thus induce disorder through ge- 
ometry. 

The resulting measurements show clear yield stress be- 
havior and can be excellently fit by the Herschel-Bulkley 
model, yielding for all bubble sizes (3 w 0.4, which is 
markedly lower than the 0.70 found for the drag force in 
ordered lanes above, and close to the 0.36 extracted from 
the velocity profiles (see Fig. [7^-b). The observed stress 
plateau at low strain rates increases with increasing bub- 
ble radius, contrary to the intuition that the yield stress 
is set by the Laplace pressure and should hence scale in- 
versely to the bubble radius. We tentatively attribute 
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this to the deformation of the bubbles through capillary 
effects, which are larger for larger bubbles and hence lead 
to a relatively larger contact size between the bubbles. 

In order to further establish a connection between the 
rhcometrical data and the model, we now turn to a ge- 
ometry with a large inner wheel to increase the measured 
signal (r, = 5 cm and r a = 7 cm), and measure the torque 
exerted on the inner wheel by a bidisperse foam with 
the same bubble sizes as in the linear shear experiment. 
We obtain a clear confirmation that indeed the disorder 
changes the power law scaling of Fbb'- we again repro- 

ducibly measure Herschcl-Bulklcy behavior with power 
law index /3 ss 0.40, as can be seen in Figs. [7]jc). 

To convert torques to F bb , we divide the torque by the 

number of bubbles and . Since our outer rough bound- 
ary forces the bubble velocity to zero, we can rescale the 
angular frequency to the dimcnsionlcss velocity difference 
T]Av/cr by assuming a linear velocity profile across the 
gap, decaying from tori to 0. The gap width is approxi- 
mately 9(d) and hence we can estimate Aw. We extract 
from the rheological measurements an estimate for the 
ratio k = fbw/ fbb ~ 5.5 ± 0.5. This is remarkably close 
to the value k = 3.75 ± 0.5 extracted from the veloc- 
ity profiles, given the crude estimates used in converting 
torques to bubble-bubble drag forces in the rheometrical 
data — we have oversimplified the shape of the veloc- 
ity profile in the disordered Couette rheometry, which is 
neither linear, nor rate independent. 



C. Interpretation 

The drag forces exerted on the bubbles by the top 
plate, which at first sight might be seen as obscuring the 
bulk rheology of the foam, enable us to back out the ef- 
fective inter-bubble drag forces and constitutive relation 
of foams from the average velocity profiles. To further 
appreciate this fact, note that our model yields linear ve- 
locity profiles regardless of the exponent (3 if the body 
force due to the wall drag is zero. This is consistent with 
earlier measurements by Wang et al. [IS], where essen- 
tially linear flow profiles were found for bubble rafts, i.e., 
in absence of a top plate. 

By comparing the results obtained from the velocity 
profiles with the rheometrical measurements, we note a 
remarkable difference between the scaling of the bubble- 
bubble drag forces at the bubble level, which we have 
mimicked by strictly ordered bubble rheology, and the 
scaling of the averaged forces at the bulk level, which we 
have extracted from the velocity profiles and confirmed 
by rheometry: we find Fbb ~ (Av) ' 70 at the bubble level 
and F bb ~ (Au) 36 at the bulk level. 

We speculate that this is closely connected to the non- 
affine behavior of the bubbles [51 [331 03]: close to the 
jamming transition, the effective viscosity of the foam 
becomes anomalously large due to the fact that bubble 
motion is much more complicated than if the bubble mo- 
tion would have been affine. i.e., where the bubbles follow 



10 ! 
10' 
10° 
£10"' 
°-10 J 
10" 
10" 



1.0 



X 0-5 



0.0 



(a) 


■ 41 = 016 a 
ft 4 at= 23.15 s 

A 


d 

t ■>? 




6-4-2 


2 4 f 


(c) 




At = 0.46s 


™ * 5-0.6 

• 5 -07 

♦ E-0J 


0.0 


0.5 1.0 



10' 
10' 
10 s 
10' 

10 ■' 

10 s 
10 J 



(b) 



1 ■ At = Q.AB s 
- At = 23,13 s 



A 



-6 -4 -2 2 4 6 



1.0 



0.5 



0.0 



(d) 




it = 23.15s . 






■ S -0.6 




• i -0,7 




A i -0.6 




+ 4 -0.9 



0.0 



<(v/v/"> 



0.5 1.0 



FIG. 8: (Color online) Dimensionless velocity distributions 
measured for W = 7 cm, no = 0.25 mm/s, and for a y-position 
19 mm away from the center of the gap. Here the averaged 
velocity equals 3.1 x 10~ 2 mm/s and the local strain rate 
equals 5.6 x 10 -3 s _1 . (a) Distribution of v x /vo for a short 
time interval (At = 0.46 s, black squares) and longer time 
interval (At — 23.15 s, red circles). For the averages of the 
dimensionless velocity distributions we find (v x /vq) « 0.125, 
independent of the time averaging interval, (b) Distribution 
of (vx/vo) 2 ^ 3 for a short time interval (At — 0.46 s, black) 
and longer time interval (At = 23.15 s, red). The aver- 
ages of the scaled dimensionless velocity distributions equal 
<(^/«o) 2/3 )o.46s » 0.205 and {(v x /v ) 2/3 ) 23 . 15s = 0.245. The 
significance of this is that (v x /vo) 2 ^ 3 ~ 0.248, which is signif- 
icantly better approximated by the longer time average, (c) 
Comparison of {v x /vq)^ and (v x /vq) 2 ^ 3 along the flow profile 
for At = 0.46 s, and for four values of £ as indicated. The 
best linear relation is obtained for £ w 0.80. Dotted vertical 
line indicates the averages shown in panel (b). (d) Same as 
(c), now for At = 23.15 s. The best linear relation is obtained 
for £ w 0.72. 



the imposed shear |35j . This picture is corroborated by 
recent simulations on the bubble model [5] , where one re- 
covers this "renormalization" of the drag force exponent 
[401142 , 43J. The precise microscopic mechanism, though, 
is far from understood. 

One may wonder why the modification of the exponent 
of the drag force law is strong for the inter-bubble forces 
but weak or essentially absent for the bubble wall drag 
forces. We have no definite answer, although we are fairly 
confident that the bubble-wall drag forces indeed are not 
modified. We base this assertion on explorations of the 
bubble trajectories, described below. 

If we assume the Bretherton expression, Eq. Q, to 
be the correct expression that gives the instantaneous 
bubble-wall drag force as a function of the instantaneous 
bubble velocity, our claim is that the averaged bubble- 
wall drag forces scale similar to the individual bubble- 
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wall drag force: 



((v/\v\) x \v\ 2/3 ) « (v x ) 2/3 



(15) 



Hence we claim that the time averaged bubble wall drag 
force is proportional to (v x ) 2 ^ 3 , which is the expression 
we employ in our model to estimate Fb w . In other words, 
we can interchange the order of taking time averages and 
"raising to the power 2/3". 

To check this, we have performed accurate bub- 
ble tracking and calculated and compared (Fb x ) = 
{(v/\v\)M 2/3 ) and K) 2/3 [13 EES- In Fig.[8k-b we show 



examples of distributions of both (Fbx) and 



V2/3 



, based 



on short and long time velocity estimates at a fixed posi- 
tion in the cell. For long times these distributions are nar- 
rower and have less weight around zero. For the examples 
shown in Fig. [8^-b, the averages of the dimensionless ve- 
locity distributions equal (v x /vq) w 0.125, independent of 
the time averaging interval. Hence, (v x /vq) 2 ^ 3 w 0.248. 
The averages of the distributions of (v x /vq) 2 / 3 , taken 
over different time intervals, depend now on this time 
interval and approximate (v x /vq) 2 / 3 ~ 0.248 better the 
longer the time interval is: we find ((v x /v ) 2 / 3 ) A6s ~ 
0.205 while ((v x /v ) 2 / 3 ) 23 .i5 S = 0.245. Since the drag 
force model deals with (long) time averages, the improve- 
ment of the agreement with time is encouraging. 

The connection between (Ff, x ) and (i> x ) 2 ' 3 can be 
probed in more detail by plotting (Ff, x ) as function (v x )^ 
for a range of strain rates, and estimating for which value 
of £ these two quantities are proportional. The data in 
Fig. [8}; shows that for short times, a value of £ ~ 0.80, 
significantly different from 2/3, leads to the best corre- 
lation, while for longer times — |8jd) — , the best value 
is £ ss 0.72. Therefore, the longer the time interval, the 
closer £ approaches 2/3. The underlying reason is that 
for increasing time intervals, the distribution of v x /vq 
becomes narrower and narrower and peaked away from 
zero, and thus we indeed can interchange the order of 
taking time averages and "raising to the power 2/3". 

The bubble-bubble drag forces, on the other hand, 
involve velocity differences, and even at long times we 
expect their probability distribution to have significant 
weight around Av = 0. The situation is then qualita- 
tively similar to that shown in Fig.[8]for short times, and 
a change from local to global exponent appears reason- 
able. Unfortunately, testing this explicitly in our data 
for the bubble trajectories has proven to be prohibitively 
difficult, not only because velocity differences are smaller 
and more noisy than velocities, but also since bubble 
contacts are very hard to establish unambiguously. The 
precise mechanism responsible for the "renormalization" 
that leads to the exponent /3 ~ 0.4 remains therefore 
open. 

Finally, the origin of the edge effects that prevent us 
from fitting our full experimental curves with the model 
profiles, might be due to the fluid drag near the wheels 
that was discussed in section II C. Alternatively the ori- 
gin might lie in the absence of a local flow rule near the 
driving wheels as reported in |54j . One way to resolve this 



1.0 - 



0.5 



^0.0 



-0.5 



-1.0 - 



- I 1 1 1 1 1 1 1 1 1 1 1 1 - 



2 area% defects 



r 



(a) 




G 1 



9,2 area% defects 



□STBCI5 



(b) 



_i i i i i i i i i ; i i i_ 



-12 8 



4 4 

y/<d> 



12 



FIG. 9: (Color online) (a) Velocity profiles for a monodis- 
perse, ordered foam with the crystal axis aligned with the 
wheels. Gap W = 7 cm and v = 0.083 (black), 0.26 (dark 
grey) and 0.83 (light grey) mm/s. Solid curves indicate fits 
to the model Eq. JT5| with k = 0.3, ,9 = 2/3. (a-b) Veloc- 
ity profiles for an ordered foam consisting of 2.7 mm bubbles 
for same drving velocities as main panel, to which defects are 
added in the form of an increasing area fraction of 1.8 mm 
bubbles as indicated. 



is accommodating non-local behavior in our model, for 
instance by incorporating drag terms due to next nearest 
lanes, similar to the cooperativity length introduced in 
|54j . We have not pursued this avenue. 



V. ORDERED FOAMS 

We have postulated that the disordered bubble mo- 
tion underlies the anomalous relation between the lo- 
cal bubble-bubble drag forces and the global viscous 
stresses. To corroborate this conjecture, we shear or- 
dered, monodisperse foams in the linear geometry, simi- 
lar to what was done in [18] . In this case the bubbles are 
expected to move affinely with the global shear, in which 
case one would expect the global viscous drag forces to 
scale the same as the local ones. 

We shear a monodisperse, ordered foam with bubbles 
of size 2.7 mm, produced by blowing nitrogen through 
one syringe needle at fixed flow rate, at a gap IF of 7 cm 
at «o = 0.083,0.26 and 0.83 mm/s. We recover the rate 
independent and strongly shear banded velocity profiles 
reported in [TS] (see Fig. [9]). As in the case of the bidis- 
perse foams, we fit model profiles to our experimental 
data. For our model to yield rate independent velocity 
profiles, the drag forces need to balance in the same ra- 
tio for all driving velocities. This can only be achieved 
if /3 = 2/3 since we have already confirmed with rheome- 
try that the exponent governing bubble- wall drag is 2/3. 
Indeed we find that the experimental profiles are best fit 
by model profiles if one fixes k — 0.3 and j3 — 0.67 ± 0.05 
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|, see Fig. [9] 



Disorder 



In our experiment, the complex bubble motion is 
closely connected to the anomalous scaling of the bubble- 
bubble drag force, which in turn is reflected in the ob- 
served rate dependence of the velocity profiles. We can 
thus investigate for which levels of disorder the rate de- 
pendence of the velocity profiles occurs by gradually in- 
creasing the disorder, starting from a monodisperse foam. 

To this end we record velocity profiles in a monodis- 
perse foam made of 2.7 mm size bubbles in which we 
gradually increase the area fraction of smaller (1.8 mm) 
bubbles. After mixing the two species we measure ve- 
locity profiles at v = 0.083,0.26 and 0.83 mm/s. We 
already observe the occurrence of rate dependent veloc- 
ity profiles for small quantities of defects, sec inset (a) 
and (b) of Fig. [9j and by visual inspection, we already 
see the swirling patterns, typical of our 41/59 bidisperse 
foam, occurring at 2 % disorder. These findings indicate 
that rate independent flows are in fact limited to a nar- 
row region close to the almost singular case of completely 
ordered foams. 



VI. ROLE OF THE PACKING FRACTION 

In this section we will discuss linear shear experiments 
where we will vary both the packing fraction (or wetness) 
of our foam <p as well as the applied strain rate, to inves- 
tigate the flow behavior of these foams as a function of 
density. In particular, we will closely approach the jam- 
ming transition, located at 4> c « 0.84. This allows us to 
test our drag force balance model over a wide range of ex- 
perimental situations. Our main findings are that, first, 
the scaling exponent (3 appears to be independent of 
and second, that the pre-factor k is our model (Eq 
varies as 1/(0 — (f> c ) where 4> c ~ 0.84. 
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A. Varying and measuring 4> 

In order to vary 0, we vary the vertical gap between 
the glass plates and the bulk solution between 3 and 0.2 
mm. We do this by adding or retracting fluid from the 
reservoir. For large gaps the bubbles get stretched in the 
vertical direction, and share large deformed facets — the 
foam effectively becomes dry. For small gaps the bubbles 
acquire a pancake- like shape, close to purely disc like 
in the horizontal plane, with only small facets between 
neighboring bubbles — the foam effectively becomes wet. 

To create a homogeneous gap between the liquid sur- 
face and the glass plate, we place additional supports 
under the glass plate to prevent sagging of the top plate 
during the runs. We monitor the gap width with a Mi- 




FIG. 10: Image manipulations leading to a definition of <j>. 
Left: Raw image. Center: Raw image with reconstructed 
bubble areas superposed. Note the good agreement. Right: 
Final binarized image from which packing fraction is deduced. 



tutoyo digital depth gauge. If the gap becomes smaller 
than 0.2 mm the bubbles unjam |56l 157] . 

We find that in the linear shear cell the accessible range 
in <fi is 0.86 < <j> < 0.97. If we stay between these limits 
the system we study is jammed and quasi two-dimensi- 
onal. It should be noted that for the runs performed at 
fixed wetness, discussed in the previous sections, we find 
4> = 0.965 ±0.005, in reasonable agreement with previous 
reports on the maximum <fi that can be obtained in our 
type of setup [55] . 

Determining the liquid fraction is not trivial, since var- 
ious horizontal cuts through the bubble layer will yield 
different values 59J . We choose our lighting of the bub- 
bles such that the contacts between adjacent bubbles arc 
optimally resolved. We then extract <f> through image 
analysis, as illustrated in Fig. [10] We first binarize the 
images, after which both the bubble centers and the in- 
terstices appear bright. We remove the interstices by 
morphological operations. We then invert the binarized 
image and fill up the remaining bubble contours. We have 
checked that the resulting bright disc optimally matches 
the original bubble contour, see Fig. 



10 We then cal- 



culate the ratio of white pixels over the total number of 
pixels and hence obtain a reasonable estimate of 4>. 

Now that we have obtained good estimates of the pack- 
ing fraction <j), we can probe the role of the wetness in 
setting the flow. We first, in section |VIB[ briefly dis- 
cuss a local probe of the non-affine motion, which shows 
that the bubble motion becomes increasingly non-affine 
when the wetness is increased. We then investigate the 
variation of the flow behavior with (j), using our model 
Eq. (10 1. We first establish, in section VIC that the 



exponent (3 does not vary with cj> — surprising, give the 
varying degree of non-affinity. We then find, in section 
|VID| that the force pre-factor k varies strongly with <f> 
and vanishes at 4> c ~ 0.84 as 1/(0 — 0c)- 
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FIG. 11: (Color online) Displacement angle distributions 
P(a) for runs for which vo = 0.26 mm/s, W = 5 cm 
and At — 0.46 s averaged over the shear banded region ( 
< y < W/3 and 2W/3 < y < W) for the range of packing 
fractions as indicated. 



B. A local measure of the non-affine bubble 
motion: P(a) 

A crucial feature of deformations found in simula- 
tions of packings of frictionless discs near jamming is the 
strongly non-affine nature of the particle (bubble) motion 
[5J |35J [371 BI] • Recently, a simple local probe of this affin- 
ity was introduced by Ellenbroek et al. who performed 
simulations of soft frictionless discs [44 . Defining the 
displacements of contacting particles i and j as iij and 
Uj, and the vector that connects the centers of particles 
i and j as fj , the relative displacement angle a was de- 
fined as the angle between fj and Ui — Uj. In other words, 
a = 0° corresponds to particles moving away from each 
other, a = 180° corresponds to particles moving closer, 
and a = 90° corresponds to particles sliding past one 
another. 

The probability distribution P(a) was found, for shear 
deformations in particular, to be well fitted by a (period- 
ically extended) Lorentzian peaked around 90° [4"4l ISO]. 
The width of the peak scales with distance to jamming 
— at jamming, P(a) approaches a delta function peaked 
at a = 90°. 

Of course, in our experiment we have flow, and we can- 
not determine deformations in linear response. Moreover, 
our system is not homogeneous. Nevertheless, as a coarse 
measure of the degree of non-affine bubble motion, which 
we claim underlies the anomalous scaling exponent (3 in 
disordered systems, we have calculated P(a) focussing on 
finite time displacement fields (vq — 0.216 mm/s, W = 5 
cm, At = 0.46 s). 

In Fig. [IT] we show P(a) averaged over the regions 
< y < W/3 and 2W/3 < y < W where most of the 
flow takes place, and averaged in the x-direction over 50 
mm in the center of the cell. We limit ourselves to this 
region, because, in particular for the wet runs, there is 
hardly any flow in the center region of the cell and the 




FIG. 12: (a) variance in k values for all six runs performed at 
4> = 0.905 (grey squares) and 4> — 0.925 (light grey squares). 
The variance at <f> = 0.965 (black squares) is data from 
Fig.gf). A clear minimum can be observed around /3 = 0.38. 



peaks in P(a) are less pronounced in this region. We find 
that, analogous to what is found in simulations [44] . the 
distributions become increasingly peaked around a = 90° 
for increasing wetness. Moreover, the distributions are 
well fit by the same Lorentzian fit that also captures the 
numerical displacement fields well [00] , 

Hence, this simple measure of non-affine motion 
strongly indicates that the degree of non-affinity increases 
for wetter foams. We believe that this is the first exper- 
imental measurement of this distribution that shows the 
proximity of the jamming transition. Detailed studies of 
the role of the local strain rate or the time interval over 
which displacements are measured arc deferred to later 
work. 



C. Variation of the exponent j3 with <j> 

We now investigate the validity of applying the drag 
force balance model with a fixed (3 — 0.36 for varying </>. 
The microscopic exponent 2/3 which governs the flow of a 
bubble past a wall appears to be independent of the par- 
ticularities of the foam flow [3"0l IBTj . On the other hand, 
it is not at all obvious that /?, which governs the averaged 
bubble-bubble drag forces, does not depend on <fi. As we 
have seen, (3 is set by the disorder in the system and the 
non-affine bubble motion that occurs in conjunction with 
that, and as we have shown in the previous section, the 
degree of non-affinity varies substantially with <\>. 

To see if (3 indeed depends on the foam density we per- 
form two additional scans over the same six shear rates 
as employed in section III for a bidisperse foam at a gap 
width W — 7 cm, while first fixing <f> = 0.905 ±0.005 and 
then 4> — 0.925 ± 0.005. We look for a minimum of the 
variance in k over the six velocity profiles as a function 
of (3 (see grey and light grey squares in Fig. 12 1. We ob- 
serve that the model fits best to all six runs performed 
at <f> = 0.905 for a = 2/3, /3 = 0.38 ± 0.05 (see Fig. [l2l 
and k — 7.5, whereas the model best matches the runs 
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FIG. 13: (Color online) Velocity profiles in linearly sheared 
foam at fixed driving rate (vo = 0.26 mm/s), for <f> varying 
between 0.855 < <j> < 0.975 at W = 7 cm. Fits are solutions 
to linear drag force balance model with a = 0.67 and /3 = 0.36 
fixed, k is extracted from the fits and plotted in Fig. [14] as a 
function of d> — d> c . 



performed at = 0.925 for a = 2/3,(3 = 0.39 ± 0.05 
(see Fig. 12 1 and k = 5.8, thus strongly indicating that 



within our range of accessible liquid fractions (3 seems to 
be constant. For comparison, we include the variance for 
the runs described in section III B that was plotted in 
Fig. |4jf). Remarkably, (3 remains a constant with vary- 
ing <f> while the degree of non-affinity varies. While we 
do not pretend to understand this, we do remark that 
(3 and P(a) essentially encode different routes towards 
jamming and thus towards increasing non-affinity: (3 is 
renormalized by the increasing non-affinity as one lowers 
the strain rate 7 towards jamming, while P(a) monitors 
non-affinity as a function of density. 



D. Scaling of the force pre- factor k with (f> 

Now that we have established that (3 is independent 
of 0, we will probe the variation of k with </>. We mea- 
sure averaged velocity profiles at gap widths W = 5 cm 
and W = 7 cm and fixed vq = 0.26 mm/s (the third 
slowest driving velocity), for packing fractions varying 
between cf> — 0.855 and tfi — 0.975. The velocity profiles 
for W = 5 cm are plotted in Fig. [13J and are seen to be- 
come increasingly shear banded as we approach (f> c |64| . 
This trend is reflected in the increase of k as we approach 
4> c . We obtain k by fitting solutions of our drag force bal- 
ance model with a — 0.67 and (3 = 0.36 to these profiles. 



The resulting fits are shown as red lines in Fig. 13 and 
fit the data well. 

<t>r, With <t> r 
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FIG. 14: (Color online) Scaling of k with A^ = (f> — (f> c - Trian- 
gles: data obtained from fits depicted in Fig. |13| where W = 5 
cm. Squares: data for gap of 7 cm. Large squares correspond 
to runs at vo =0.26 mm/s from Fig. 12 Solid line: 0.45/A^. 
Inset: same data on log-log scale. 



the theoretically predicted and experimentally measured 
value of the unjamming packing fraction: <j) c = 0.842 
[40, 62, 63 . In good approximation we obtain that 



k cx !/(</> - 4> c ) 



(16) 



We can tentatively explain the observed scaling of k 
with a simple argument based on the sizes of the facets in 
the foam. At fixed <fr, our drag force balance model yields 
a value of k that sets the relative influence of the bubble- 
wall drag with respect to the bubble-bubble drag and 
which we have conjectured to be given by k oc fbw/fbb- 
As we have already discussed, fb w oc ar c with r c the 
radius of the flattened contact between the bubble and 
the wall and fbb oc crn^., with n c the radius of the flattened 
contact between neighboring bubbles. Thus we expect: 



k cx r c j ' n 2 c 



(17) 



While r c is set by the buoyancy and hence does not vary 
strongly with the gap distance between glass plate and 
liquid surface — only becoming slightly smaller as the 
bubbles get stretched at large gaps — n c is strongly de- 
pendent on the gap size and hence on the packing fraction 
of the foam. 

The size of k c should depend on the deformation (also 
called the overlap) S£ as |25] : 



1/2 



(18) 



Similar to simulations of two-dimensional frictionless 
discs [3H1 H3 we can relate the overlap to the packing 
fraction <h: 



5£ cx A(j). 



(19) 



In Fig. 14 we plot A: as a function of 



Simple substitution of this result into Eq. (17 1 yields 
fccxr c /^«l /8£ = 1 /(A0) , (20) 
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which is fully consistent with our experimental results, 



see the solid line in Fig. 12 



Note that in the above we have only focussed on 
the radius of the deformed facets. A proper anal- 
ysis would include the size of the Plateau border 
around the contact, which is where the dissipation also 
occurs [210 EH]- For instance, in [ST] the bubble- 



wall drag force scales as F 



OC 



Ca° 



64, 



-0.26 



and 



a proper treatment would entail such analysis, even 
though the functional dependence on the Plateau border 
size is always weak. Moreover, in all of these works, the 
functional dependence of the drag force with cj> is smooth 
around <p c and hence will not influence the observed scal- 
ing around that point. 



VII. DISCUSSION AND CONCLUSION 

We have measured velocity profiles in linearly sheared 
quasi-two-dimensional foams in the liquid-glass configu- 
ration. We find that bidisperse, disordered foams ex- 
hibit strongly rate dependent and inhomogeneous (shear 
banded) velocity profiles, while monodisperse, ordered 
foams are also shear banded, but essentially rate inde- 
pendent. We capture these findings in a simple model 
that balances the drag forces in our system. The scaling 
forms for these drag forces are verified by independent 
rheological measurements. Finally, we apply our model 
to velocity profiles obtained for foams at varying pack- 
ing fraction, and measure and describe the scaling of the 
inverse foam consistency with packing fraction. 

This work raises several questions. First, can the 



difference between the local bubble-bubble drag force 
scaling and the global (averaged) bubble-bubble drag 
force scaling be understood theoretically? This differ- 
ence in scaling exponents appears similar to the change 
from local drag forces to global rheological laws, ob- 
served in simulations of (variants) of the bubble model 
[3 SOI ED E31 [5S], but a precise connection is lacking 
at present. Closely connected, is our scenario an exam- 
ple of a general route by which aspects of the ubiqui- 
tous Herschcl-Bulkley (power law) rheology observed for 
a wide range of disordered materials can be rationalized? 

Second, how robust are our experimental results? For 
example, would similar flows in Hele-Shaw cells behave 
differently, as suggested by the results of Debregeas [15] ? 
We also wonder if our model is able to capture shear 
banded flows in Couette geometries, where the curvature 
plays an important role, in particular since the foam has 
a finite flow threshold [66]. Third, can similar phenom- 
ena and models as described here be extended to three 
dimensional flows of foams and emulsions — where flows 
in the latter can be captured by confocal imaging and 
MRI [271 [671 EE]? Fourth, how should our local mod- 
els be compared to the non-local effects recently discussed 
for emulsion flows [531 167j ? 
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